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ABSTRACT 

New high quality Keck and ESO images of PKS 1830-211 are presented. Applying 
a powerful new deconvolution algorithm to these optical and infrared data, both images 
of the flat spectrum core of the radio source have been identified. An extended source 
is also detected in the optical images, consistent with the expected location of the 
lensing galaxy. The source counterparts are very red at I — K ~ 7, suggesting strong 
Galactic absorption with additional absorption by the lensing galaxy at z = 0.885, and 
consistent with the detection of high redshift molecules in the lens. 
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Subject headings: cosmology: observations, gravitational lensing — quasars: individual 
(PKS 1830-211) — infrared: general — techniques: image processing — methods: 
data analysis 

1. Introduction 

The bright radio source PKS 1830-211 (Subrahmanyan et al. 1990; hereafter S90, Jauncey 
et al. 1991) has attracted much attention as the most detailed example of a lensed radio ring. 
Among the classically-lensed QSOs, its short time delay of 44 days (van Ommen et al. 1995) 
and clean lens geometry make it a good candidate for measuring Ho. The lens, a gas rich galaxy 
at z=0.89, was discovered in the millimeter via molecular absorption (Wiklind & Combes 1996), 
which is seen towards only one of the two flat spectrum hot spots (Wiklind & Combes 1996, Frye 
et al. 1997). The presence of a nearby M-star as well as heavy extinction along the line of sight 
(b=-5.7 degrees) had until now hampered the identification of the lens and the source. In this 
paper we describe how the MCS deconvolution algorithm (Magain, Courbin, & Sohy 1998) was 
used to detect the counterparts of this bright radio ring in deep Keck optical and infrared images. 

2. Observations - Reductions 

Near IR J (A c =1.25 micron) and K' (A c =2.15 micron) images were obtained on the nights 
of 1997 April 14 and 15 with the IRAC2b camera on the ESO/MPI 2.2m telescope, which 
uses a NICMOS3 256x256 HgCdTe array. The good seeing (0"6-0'.'7) and the fact that a good 
sampling is desirable for deconvolution, led us to choose the smallest pixel size available on this 
instrument, i.e., 07151, resulting in a total field of view of 38" on a side. The data were taken 
and reduced exactly in the same way as in Courbin, Lidman & Magain (1998). Several standard 
stars were observed during the night. The standard deviation in the zero points was less than 0.02 
magnitudes for both J and K' . The IR magnitudes reported in this paper are on the JHK system 
of Bessell & Brett (1988). 

Near-IR Keck I data were obtained on the night of 1994 April 5 with NIRC (Matthews 
& Soifer, 1994). Five 80 second K images were combined using the sigma clipping algorithm 
available in IRAF. The pixel size is 07157, similar to that used with IRAC2b. The images were 
obtained under marginally non-photometric conditions, with a zero point uncertain by about 
0.1 magnitude. Due to the crowding of the field, and the low number of frames available, sky 
subtraction was imperfect but did not influence dramatically the quality of the data analysis. 

Six dithered /-band images were obtained during the full moon night of 1997 June 15 using 
the Keck n telescope and LRIS (Oke et al. 1995). The CCD detector was a Tektronix 2048x2048 
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with a pixel size of Of! 215. The individual exposure times were restricted to 3 minutes in order 
to avoid saturation of the brightest stars in this extremely crowded field. The background was 
very high. The images were bias subtracted and flat fielded in the standard way. All of them have 
a seeing close to Of! 8. No standard star was observed, but a flux calibration could be performed 
relative to an /-band image taken at La Silla with the 0.9m telescope on the night of 1997 April 
15. 



3. Image deconvolution 

The MCS deconvolution code (Magain et al. 1997) was applied to all the images. Due to the 
low signal-to-noise (S/N) and the numerous bad pixels in single IR images, these were medianed 
in groups of 9 dithered and sky subtracted frames. The resulting images have better S/N and 
cosmetics. Two infrared-bright nearby stars, although outside the field of view, produced scattered 
light in the field, forcing us to reject a fair fraction (~ 40 percent) of our observations in the K' 
band. One of the culprits was the known source IRAS 18306-2106. Two stacks were obtained in J 
(total exposure time of 1920 sec) and four in K' (total exposure time of 1440 sec). Only one such 
stack was obtained for the IR Keck images since we had only 5 dithered frames to combine (total 
exposure time of 400 sec). 

3.1. Application of the MCS code to the present data 

The deconvolution process is the same as in Courbin et al. (1998). We chose a pixel scale in 
the deconvolved images that is a factor of two smaller than the pixel scale in the original data, to 
insure that the sampling theorem is satisfied in the deconvolved images. After deconvolution, the 
resolution is fixed to 2 pixels Full- Width-Half-Maximum (FWHM) for all data. The corresponding 
final pixel scales and resolutions on the plane of the sky are given in Fig. 1. 

We constructed the infrared Point-Spread-Function (PSF) from the relatively isolated, bright 
star labelled H in Djorgovski et al. (1992; hereafter D92). In the optical images, the stars 
mentioned in D92 are all saturated. Consequently, the PSF was constructed from 4 fainter stars 
about 30" away from the radio ring. Crowding in the optical field made the quality of the PSF 
rather poor, especially in the outer wings of the hexagonal Keck PSF, introducing systematic 
residuals into the deconvolved image and residual maps (see section 3.2). 

In each band, all the frames available were deconvolved simultaneously, i.e. 6 in /, 2 in J, 4 in 
K' and 1 in K. In other words, the output of the procedure is a unique deconvolved image which 
is simultaneously compatible with all the images included in the data set. The deconvolved image 
is given as a sum of point sources with known Gaussian shape and a diffuse background smoothed 
on the length scale of the final resolution chosen by the user. The photometry and the astrometry 
of the point sources are also obtained as byproducts of the deconvolution and are provided in 
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Tables 1 and 2. 

3.2. Quality check of the deconvolution 

Many deconvolution algorithms generate the so-called "deconvolution artefacts" as well as 
noise enhancement. Even if the MCS code does not produce artefacts, it still has to accomplish 
the difficult task of deconvolving the image from imperfect and noisy data. An objective criterion 
has therefore been established to check the quality of the results. A natural product of the 
MCS code is the RM (Residual Map), which is the difference between the deconvolved and real 
images divided by the standard deviation of each pixel. Whenever the simultaneous deconvolution 
capabality of the method is used, as in the present case, the program returns one RM for each 
frame in the data set. A perfect deconvolved image should match the data at best in the sense of 
the x 2 an d leave flat RM with a mean value of 1 (one standard deviation) all over the field. 

RMs guide the user in his/her choice of the smoothing to be applied to the image in order to 
avoid local under- or over-fitting of the data. Still, on the basis of the RMs, the user can constrain 
the number of point sources to be involved in the deconvolution. For example, missing a bright 
point source (i.e., with high S/N) will result in a "hot spot" well above the critical value of 1 in 
the RMs. This is no longer true for very faint sources. In this case, no discrimination can be made 
between an extended object and a point source, so the faint point sources are modelled as part of 
the deconvolved background. We therefore always choose the minimum number of point sources 
needed in order to produce acceptable RMs. Thus, different images of the same field can have 
different aspects, depending on whether individual objects were modeled as point sources or not. 

Objects near the frame edges are not well-fitted, especially in crowded fields such as these 
(e.g., objects labelled 1 to 4 in Fig. 1). Object 1 leaves particularly significant structures in the 
RMs of the /-band data, either when modeled as a point source or as part of the background in 
the deconvolved image. We therefore conclude that it is extended or that it is a very strong blend 
of point sources. 

In order to run the MCS algorithm, the user has to provide initial positions and intensities for 
all points sources. The astrometric la error bars were estimated as the dispersion of the results 
of several deconvolutions using different initial conditions. The typical astrometric accuracy is 
(X'05 for the brightest point sources, while it can be as much as (y/2 for the faintest ones. Since 
the peak intensity of a point source is allowed to be different in each data frame, the photometric 
error bar is simply the dispersion of the peak intensities found after simultaneous deconvolution 
of a whole data set. In order to quantify the additional errors introduced by a poor knowledge of 
the PSF in the Keck /-band, the deconvolution program was run using 5 different PSFs computed 
from different stars in the field. 



- 5 - 



4. Results and discussion 

In Fig. 1 and Plate 1 we present both the raw and deconvolved images in /, K' and K, 
with a resolution of the order of 0?15-0?20. A red point source is clearly detected at the position 
expected for the NE radio source of PKS 1830-211. Another red object is observed close to the 
position of the SW radio source of the lensed system, but the extended nature of the source and 
the poor quality of the PSF do not allow to sort out its detailed morphology. The photometry and 
astrometry of the field are presented in Tables 1 and 2 and in Fig. 2, along with our estimates of 
the la error bars. 

The red star in D92 dominates the total flux in the /-band, while in the infrared the 
component near the NE radio source takes over. With our high S/N we can show that its shape is 
compatible with a point source, and its color, I — K = 6.9, is much redder than any "normal" star 
(e.g., Koornneef, 1983). The /, K', and K positions are all within the la radio error bars. In Fig. 
2 all objects are aligned with respect to the M-star so that we can compare the different positions 
measured for the NE and SW components at optical and IR wavelengths. 

Molecular absorption towards the two lensed images was spatially resolved at millimeter 
wavelengths, and the separation between the lensed images was found to be 0"98 (Frye et al. 1997). 
In our optical and IR images the SW component is 0.61 ± 0.13" and 0.85 ± 0.09" respectively, 
away from the NE component. A plausible explanation for the apparent positional shift between 
the optical, IR and radio positions is that the SW component is a blend of two objects: the 
lensing galaxy and the heavily reddened SW component seen in the radio images. This has been 
very recently confirmed by HST NICMOS3 observations by the CASTLES group (Kochanek et 
al. 1998) which show the lensing galaxy in the center of the "hole" of the radio ring, as predicted 
by LensClean modeling of the system by Kochanek & Narayan (1992). Note also the very good 
general agreement between the HST images and our deconvolutions, although the latter were 
carried out without prior knowledge of the HST data. 

The flux ratio between the two lensed images of the source is 0.973 ± 0.016 at radio 
wavelengths (S90), and about 18.2 in K (3.15±0.2 magnitudes). The combination of a reddened 
SW radio source plus blue lens can explain the large flux ratios. Both NE and SW components are 
reddened, making PKS 1830-211 another good example of a dusty lens, along with MG 0414+0534 
(e.g., Annis & Lupino 1993) and MG1131+0456 (Larkin et al. 1994), the mean galactic extinction 
being far below the values obtained for PKS 1830-211. At the lens redshift, the observed i-T-band 
corresponds to a central wavelength of 1.17 [i in the restframe. For the galactic reddening curve, 
the extinction in this band is ^.(1.2^) = 0-37 Ay (Sandage & Mathis, 1981). Thus, the differential 
extinction is about 8.5 magnitudes in the F-band and implies E(B — V) = 2.75 ± 0.2, in good 
agreement with the value of 2.4 independently derived from the HST observations (Kochanek et 
al. 1998). 

If the total extinction is comparable to the differential extinction between the two images, 
then the source is attenuated in the observed if-band by about 3 magnitudes. On the other 
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hand, considering simple SIS or point-mass models, the magnification by the lens may be similar, 
or somewhat less, maybe 2 magnitudes. So, a reasonable value for the unobscured un-lensed K 
magnitude of the source may be around 14. A typical quasar has V — K ~ 3. So the unobscured 
V magnitude of the source may be about 17, which is bright, but not unreasonable for z ~ 1 — 2. 

5. Conclusion 

The main result of our study is the detection of the optical and near-IR counterpart to the 
NE radio source of PKS 1830-211 and the possible detection of the SW component and lensing 
galaxy. However, our SW component candidate might be the lensing galaxy alone or, given the 
crowding in the field, a red galactic object almost coincident with the position of the SW radio 
source. The hypothesis of a demagnified third image of the source (S90) between the 2 main 
lensed images is unlikely as in such a case extinction of the lens would have made it visible in the 
IR. Furthermore, the IR centroid of the SW component would have been shifted towards "object 
E" rather than towards the radio position of the SW component. 

The higher contrast between the NE component and nearby M-star in the infrared make 
near-IR spectroscopy necessary for finding the source redshift. Deep, high resolution near-IR 
imaging is needed to reveal the exact nature of the faint SW component. However, even at 
the highest resolution attainable, (X'2- (f.'3 in the IR with HST (in particular in K where the 
SW component of PKS 1830-211 is best visible), deconvolution will be essential to discriminate 
between the SW component candidate, the lensing galaxy and additional faint galactic stars. 
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We also thank Alfred Rosenberg who provided the ESO /-band frame used for the photometric 
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Table 1: Summary of the photometry. 





M-star 


NE Comp. 


SW Comp.±Lens 


Limiting Mag. 


I(Keck) 


19.3 ±0.1 


22.0 ±0.2 


22.3 ±0.3 


24.0 


J(ESO) 


18.1 ±0.1 


18.7 ±0.3 


> 20.5 


20.5 


K'(ESO) 


17.3 ±0.3 


15.8 ±0.2 


19.0 ±0.4 


19.1 


K(Keck) 


16.6 ±0.2 


15.1 ±0.1 


18.2 ±0.2 


21.3 



Note. — Results are given for the M-Star, the NE component, and the SW component of the lensed source plus 
lensing galaxy. The limiting magnitude (point sources) for each of the bands is given in the final column. 



Table 2: Summary of the astrometry. 





NE Comp. 


SW Comp.+Lens 


x(/) 


±0.08 ±0.1 


+0.48 + 0.1 


y(i) 


-0.70 ±0.1 


-1.15 + 0.1 


x(J) 


-0.07 ±0.1 




y(J) 


±0.54 ±0.1 




x(K') 


+0.06 ±0.07 


+0.54 + 0.07 


y(K') 


-0.52 ±0.07 


-1.16 + 0.07 




+0.06 ±0.05 


+0.59 + 0.05 


y{K) 


-0.54 ±0.05 


-1.20 + 0.05 



Note. — Results for the NE and SW components of the lensed source plus lensing galaxy are given in arcseconds, 
relative to the M-star, together with their la error bars, 




Fig. 1. — Four arcseconds field around PKS 1830-211 observed in the I, K' and K bands at Keck 
and ESO observatories. [Top left] Stack of 6 /-band frames with a pixel size of (X'215 and seeing of 
(y/8. [Top middle] Mean of 4 temporary stacks (see text) obtained in the K' band with the ESO/MPI 
2.2m telescope. The pixel size is (X'151 and the seeing is G/.'7. [Top right] Mean of 5 K-baxi& images 
obtained with Keck I and NIRC. The pixel size is 0'.'157 and the seeing is 0?7. [Bottom right] 
Simultaneous deconvolution of the 6 /-band frames: resolution of 07215 and pixel size of 071075. 
[Bottom middle] Simultaneous deconvolution of the 4 K' images: final resolution of 07151 and pixel 
size of 07075. [Bottom right] Deconvolution of the mean of 5 NIRC images: resolution of (y/157 
and pixel size of 070785. In all the images North is up, East to the left. The M-star, the NE QSO 
component candidate and the blended SW component + lensing galaxy candidate are indicated. 
In all images the first countour corresponds to 2.5a s k y - The step between two contours is 2.5a s k y - 
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Fig. 2.— Field around PKS 1830-211 observed in the I, K' and K bands at Keck and ESO 
observatories. [Top left] Stack of 6 /-band frames with a pixel size of 07215 and seeing of 0/8. 
[Middle left] Mean of 4 temporary stacks (see text) obtained in the K' band with the ESO/MPI 
2.2m telescope. The pixel size is / .'151 and the seeing is 077. [Bottom left] Mean of 5 i^-band 
images obtained with Keck I and NIRC. The pixel size is 0/157 and the seeing is 0/7. [Top right] 
Simultaneous deconvolution of the 6 I-band frames: resolution of 0/215 and pixel size of 0/1075. 
[Middle right] Simultaneous deconvolution of the 4 K' images: final resolution of 0/151 and pixel 
size of 0/075. [Bottom right] Deconvolution of the mean of 5 NIRC images: resolution of 0/157 
and pixel size of 0/0785. In all the images North is up, East to the left. The M-star, the NE QSO 
component candidate and the blended SW component + lensing galaxy candidate are indicated. 
The horizontal line in each deconvolved image is 1" long. The cuts of the images are chosen in 
order to display the full dynamics of the image at low light levels. 
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Fig. 3. — Positions observed for the different objects detected in the optical and near-IR along with 
their la error bars, relative to the M star. The large open circles are radio positions whose radius 
corresponds to the error bars quoted by S90. The black dots show the positions measured from the 
near-IR images while the open circles show the result obtained from the /-band data. 



